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We present the analytical formulas, derived based on the halo model, to compute the cross- 
correlation between the thermal Sunyaev-Zel'dovich (SZ) eflect and the distribution ol galaxy clus- 
ters. By binning the clusters according to their redshifts and masses, this cross-correlation, the 
so-called stacked SZ signal, reveals the average SZ profile around the clusters. The stacked SZ 
signal is obtainable Irom a joint analysis of an arcminute-resolution cosmic microwave background 
(CMB) experiment and an overlapping optical survey, which allows for detection of the SZ sig- 
nals for clusters whose masses are below the individual-cluster detection threshold. We derive the 
error covariance matrix for measuring the stacked SZ signal, and then forecast for its detection 
from ongoing and forthcoming combined CMB-optical surveys. We find that, over a wide range 
of mass and redshift, the stacked SZ signal can be detected with a significant signal to noise ratio 
(total S/N ^ 10), whose value peaks for the clusters with intermediate masses and redshifts. Our 
calculation also shows that the stacking method allows for probing the clusters' SZ profiles over a 
wide range of scales, even out to projected radii as large as the virial radius, thereby providing a 
promising way to study gas physics at the outskirts of galaxy clusters. 



I. INTRODUCTION 

Galaxy clusters are the most massive, gravitationally 
bound objects in the universe, and therefore, their prop- 
erties, e.g., the cluster abundance and its redshift evo- 
lution, are very sensitive to the underlying cosmology 
including the nature of dark energy 0,0. The Sunyaev- 
Zel'dovich (SZ) effect 3 offers a unique means of finding 
clusters out to high redshift due to its redshift indepen- 
dence [see 0, for a review]. The wide-field, arcminute- 
resolution experiments of the cosmic microwave back- 
ground (CMB) such as the Planck satellite [sj, the At- 
acama Cosmology Telescope (ACT) Q, and the South 
Pole Telescope (SPT) 0, have been drawing a great at- 
tention, because such CMB surveys enable us to con- 
struct a catalog of SZ clusters under a homogeneous se- 
lection function out to high redshifts beyond z = 1. Such 
a cluster catalog can in turn be used to carry out the 
high-precision cosmological probes (sl-fiol also see refer- 
ences therein]. 

In the light of cluster cosmology, combining the SZ sur- 
vey with a multi-color optical survey offers a great syn- 
ergy, if the two cover the same region of the sky. First, 
the multi-color optical survey provides an estimate for 
the redshift of each SZ-discovered cluster from the pho- 
tometric redshifts of its member galaxies. Otherwise the 
cluster's redshift is unknown due to the redshift inde- 
pendence of the SZ signal. Second, from such an op- 
tical survey, we can construct a catalog of clusters by 
finding the concentration of gala xies in angular and pho- 
tometric redshift space (ill. Il2l]. Hence, by combining 
the SZ and optical surveys, we can obtain a more ro- 
bust catalog of clusters. Third, from the joint experi- 
ment, we can explore the cluster mass-observable scaling 



relation by combining various observables including the 
measured SZ signal strength and the optical richness (i.e. 
the number of member galaxies). We can also include the 
weak lensing measurement on an individual cluster basis 
[13I [l^ . where the lensing signal is available from the 
same optical data set, by exploring in detail the lensing 
shearing effects on background galaxy shapes due to the 
foreground cluster [see, e.g.,[l^ for a weak lensing study 
of galaxy clusters]. The use of a well-calibrated mass- 
observable relation is of critical importance to attain the 
high-precision of cluster cosmology 0, [l^ 

By cross-correlating the SZ map with the distribution 
of optically selected clusters for a sufficiently large clus- 
ter sample, one can statistically measure the SZ signals 
of these clusters - the so-called stacked SZ signal. The 
stacking method enables us to measure the SZ signals 
for halos down to relatively low masses as well as up to 
large projected radii, which otherwise are very difficult 
to measure on an individual cluster basis, because the 
signals are too small. In particular, the stacking method 
allows a tomographic reconstruction of the SZ signal as a 
function of redshift [13, El, while the CMB signal itself 
is by nature two-dimensional, reflecting the information 
projected between the last scattering surface and the ob- 
server. The stacked SZ signals are useful for studies of 
the mean cluster mass-observable relation thanks to the 
statistically-boosted signal-to-noise ratios. Furthermore, 
the stacked SZ signals can be used to study the intraclus- 
ter gas physics and to falsify and/or refine our theoretical 
understanding of galaxy clusters [l9l - [2^ . 

The initial attempts to measure the stacked SZ signals 
have been recently reported, including those by combin- 
ing the WMAP or Planck data with the SDSS maxBCG 
cluster catalog [H,!!!] and the ACT data with the SDSS 
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luminous red galaxies (LRGs) [l^], which preferentially 
reside in massive halos. Such measurements will be fur- 
ther significantly improved by combining the arcminute- 
resolution, high-sensitivity CMB surveys with overlap- 
ping wide- field, deep optical surveys, e.g., the SPT sur- 
vey with the Dark Energy Survey (DES) [l^, and the 
ACT with the Subaru Hyper Suprime-Cam (HSC) Sur- 
vey da 

Hence the purpose of this paper is to study the 
prospects of measuring the stacked SZ signals from ongo- 
ing and upcoming combined SZ and optical surveys. To 
make predictions for the stacked SZ profiles, we model 
the distribution of the hot gas contributing to the SZ 
effect d^l as well as the distribution of galaxy clusters 
according to the halo model [Slj. Our formulas should 
be useful for future stacked SZ analysis. Following [loj . 
we derive the covariances of the stacked SZ measure- 
ments for clusters binned by their redshifts and masses, 
which quantify the statistical uncertainties in the mea- 
surements. 

The outline of the paper is as follows. In Sec. [Hi we 
present our derivation for the formulas to calculate the 
cross-correlation between the SZ effect and the distribu- 
tion of galaxy clusters, and give the covariance matrix for 
its measurement. In Sec. IIIIl we forecast the prospects 
of measuring the stacked SZ signals from several ongoing 
and upcoming surveys, and discuss our results. Finally, 
we conclude in Sec. IIVI Throughout this paper, we adopt 
the best-fit flat ACDM model from WMAP 7-yr results 
[12 with flm = 0.27, rib = 0.045, h = 0.71, erg = 0.8, 
ris = 0.96. 



II. THEORETICAL CALCULATION 

In this section, following the method developed in 
and [11] [also see dl], [33 - [40j , we use the halo model to 
calculate the cross-correlation between the SZ effect and 
the distribution of galaxy clusters, and give the covari- 
ance matrix for its measurement. 



A. The SZ Effect and the SZ Power Spectrum 

The thermal SZ effect introduces the following frac- 
tional distortion to the CMB temperature measured at 
frequency v and in the angular direction of n [see e.g. 

113,111: 

AT 

— {x,n) = f{x)y{n), (1) 

where we have defined the dimensionless frequency x = 
hiy/ksT, with h the Planck constant, ks the Boltzmann 
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constant, and the frequency dependence of the effect f{x) 
is given by 

/(x) = xcoth(|) -4, (2) 

where rclativistic corrections [see e.g. l43 - l45j have been 
neglected. Note, in the Rayleigh- Jeans (RJ) limit, 
f{x) — >■ —2. The Compton y-parameter y{n) is propor- 
tional to the electron pressure Pe projected along the line 
of sight: 

y{n) = ^ [ Pede, (3) 

where ctt is the Thomson cross section, rrie is the electron 
mass, and c is the speed of light. 

Using the halo model, y{n) can be expressed as 

/ - \ o-T f cdz 1 /" , , r f 1-^ 

niec/ J H{z) (1 + z) J J 

X ^ SoiMh - M,)6%{xh -x^)Pe{\x-XhlMh,Zh), 

i 

(4) 

where H{z) is the Hubble expansion rate at redshift z; 
Sd stands for the Dirac delta function; Mi and Xi are 
the mass and center position of the i-th halo on the light 
cone, and the sum is done over all the haloes; P^Qx — 
Xh\, Mfii zji) is the pressure at position x of the electrons 
in the dark matter halo with mass M/j, redshift z^, and 
centered at Xh- Note that the 3-dimensional position 
X is specified by (Xj"-), where x is the comoving radial 
distance corresponding to redshift z on the light cone. 

By using the flat-sky approximation and the Limber 
approximation (46| . wc obtain the SZ power spectrum 

cY = cf + cf (5) 

with the one- and two- halo term contributions, Cj^'^^ 
and Cj'^'^^, given by 

where d^V/dzdVl ~ cy^ / H is the comoving volume ele- 
ment. dn/dM is the halo mass function, yi{M,z) is the 
2D Fourier transform of the Compton y-parameter profile 
of a cluster with mass M and at redshift z, and P^i^, z) 
is the linear matter power spectrum. The function W^{z) 
is defined as 

/dt} 
dM—{M,z)h{M,z)UM,z), (7) 

where 6(M, z) is the halo bias. As can be seen from 
Eq. ([6|), the SZ power spectrum arises from the contri- 
butions of halos over all the redshift range, and therefore 
contains only the projected information. 
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Below in our calculation, we would use the halo mass 
function from 1471 . specifically their Eq. (B3), and the 
halo bias from [4^, . To calculate the linear matter 
power spectrum, we use the transfer function from 50|. 
Our model for the intra-cluster gas is taken from 30j, 
where the gas is assumed to be in hydrostatic equilib- 
rium with the dark matter potential, have a constant 
polytropic equation of state, and trace the dark mat- 
ter density with a constant amplitude offset in the outer 
parts of a cluster. We use the virial overdensity [5l|, [IH 
to define the cluster mass, and truncate the gas profile at 
three times the virial radius Tvir. As implied from hydro- 
dynamic simulations, the gas profile may extend beyond 
Tvir due to various nonlinear processes such as shock heat- 
ing, feedback from star formation etc. [Sll, Hence, 
it would be interesting to study the gas properties out 
to these large radii, which can be explored by using the 
stacking analysis (see below for further discussion). 



B. The Distribution of Galaxy Clusters and Its 
Power Spectrum 

Given a cluster survey, the galaxy clusters can 
be binned according to their observed redshift Zobs 
and estimated mass Mobs- Given P(zohs\M,z) and 
P{Moi,s\M, z), i.e. the probability distribution functions 
(PDFs) of observing Zohs and Mobs for a cluster with 
true redshift z and true mass M, the probability that 
the cluster is selected to the "a"-th redshift bin with 
^obs e [^obs.min' <bs,max] and the "6"-th mass bin with 
Mobs e [M^^j, M^^^^^^J can be calculated by 



SabiM,z) 



dzQbsPizQbs\M,z) 



X / dMobsP(Mobs|M,z), (8) 

'^^obs mi„ 

where Sat is called the selection function. In the ideal 
case that both z and M can be measured accurately, 
which we would assume in the following calculation for 
simplicity, the PDFs become delta functions, and Sab 
reduces to 

Sab{M, z) = Q{z — 2;obs,min)Q(^obs,max " ^) 

X e(M - M,^b,,„,iJe(M,^b,,„,,, - M), (9) 

where Q stands for the Heaviside step function. 

With the halo model formulation, the 2-dimensional 
angular number density of galaxy clusters in the "(a6)"- 
th bin can be calculated by 



nlb{n) = / dz 



d^V 
dzdQ 



dMSab{M,: 



xY,Sd{M - M,)6l{x-x,), (10) 



where all symbols have the same meaning as before. Its 
ensemble average is given by 



-2D 

l-ab 



d^V 

dzdO. 



dMSab{M,z)^{M,z), 



(11) 



where we have used (^^5d{M — Mi)6'jj{x — Xi)) = 
dn/dM Hi. 

With similar approximations as before, the power spec- 
trum for the distribution of two cluster samples that are 
in the "(a6)"-th and "(a'6')"-th bins, respectively, can be 
calculated by 

C^Ua'^') = J dz^Prn = ^,^) W^^fc (^X-fc- (^) , 

(12) 

where W^^(z) is defined as 

W!,\iz)^^ [ dM-^iM,z)Sab{M,z)b{M,z). (13) 
Kb J dM 

As long as the redshift of clusters can be measured ac- 
curately, as we assume in this paper, only the clusters in 
the same redshift bin would correlate with one another. 



so we have Cl'^ . 



i« the 

Kronecker delta function. Note in this case, W^^fc(z) = 
if z is outside the "a"-th redshift bin. Hence, with red- 
shift information, the cluster power spectrum is free of 
projection effect. 



C. The Cross-Correlation: the Stacked SZ Profile 

With Eqs. (jH) and ([TO)) in hand, we are now ready to 
calculate the cross-correlation between the SZ signal (to 
be explicit, AT^^/Tcmb) and the cluster distribution. 
Again by adopting the flat-sky approximation and the 
Limber approximation, the cross power spectrum for the 
SZ signal and the distribution of galaxy clusters in the 
"(a6)"-th bin can be calculated by [also see[l3, for the 
method of derivation] 



W,(ah) ~ ^e,{ab) ^ ^e,{ab}^ 



(14) 



where the one- and two-halo term contributions are given 

by 



^i,{ab) 



dz^fdM^{M,z) 
dzdn J dM^ ' 



xSabiM,z)yi{M, z). 



(15) 



a 



'h,2h 
i,{ab) 



fix) 



dz^P- 
dzd^ 



k=^-,^W\{z)Wl{z). 



(16) 

This is the Fourier-transform of the so-called stacked SZ 
profile, which practically, can be obtained from observa- 
tion by averaging the SZ maps around the cluster cen- 
ters and then averaging over the orientation. Hence, the 
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center of cluster should be known a priori, e.g., from 
the position of the brightest cluster galaxy in the case 
that the cluster catalog is constructed from an optical 
survey as we have had in mind. Offsets of the "clus- 
ter centers" would cause dilution of the stacked signal 
at small-angles (high-multipoles). Since the effect can 
be straightforwardly included by using the method de- 
veloped in [l^l, we here ignore it for simplicity. We also 
notice that real clusters usually have aspherical shapes 
and substructures, which are not taken into account in 
our halo model calculation. A careful study of their ef- 
fects on the stacked SZ profile would involve comparison 
with hydrodynamical simulations (see, e.g., (H^j for their 
effects on the SZ power spectrum), which are beyond the 
scope of the current paper. 

Comparing with Eq. ([S]), one can find that, since the 
selection function Sab{z) or W^5(-j) vanishes outside the 
redshift bin of the clusters taken for the cross-correlation, 
the cross-power spectrum contains only the SZ contribu- 
tions from clusters in the redshift bin. Put another way, 
the cross-correlation method allows a tomographic recon- 
struction of the SZ signal as a function of redshift [also see 
[TtI . [isl for similar discussions] . When the redshift slice 
is sufficiently narrow, the cross-power spectrum is equiv- 
alent to the three-dimensional power spectrum through 
the correspondence of fc = ^/Xj under the assumption of 
statistical isotropy. 

Perhaps more intuitive is the Fourier-transformed 
counterpart of the cross-power spectrum, i.e., the stacked 
y-profile, defined as 
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fix) J 2tt 



Umct^a^y (17) 



where Jo{id) is the zerotli order Bessel function. The 
stacked y-profile gives the average y-profile around the 
clusters in the bin, which have similar redshift and 
similar mass, hence gives the average projected pres- 
sure profile: (/ dl Pe{r, l)){r±) = nieC^ (y) {9)/aT, where 
r± = x(-z)^'/(l + z). Under the assumption of statisti- 
cal isotropy, the average 3-dimensional pressure profile 
is considered to be spherically symmetric, and therefore 
can be recovered from the stacked y-profile through a de- 
projection process. This is analogous to the average mass 
profile of the clusters obtained from the stacked lensing 
measurement llSllSSt. 



D. 



Measurement Uncertainties: the Covariance 
Matrix 



Given a survey, the stacked SZ signal can be measured 
by constructing an estimator from the observed fields. 
Uncertainties on the estimated values can be obtained 
from their covariance matrix. Similar to the covariance 
matrix for the stacked lensing power spectrum [loj [also 
see the covariance matrix for the stacked SZ 



power spectrum is given by 



b') 



Off, 



W ^e,{ab,a'b') + ^e,{ab)^e,{a'b') 



(18) 



where /sky = fis/47r is the fraction of the survey sky, 
with ils the survey area, A£ is the width of the £ band 
used to estimate the power spectrum, and Ct with a hat 
symbol denotes the corresponding power spectrum for 
the observed fields, i.e. including noise contamination. 
Note, for simplicity, we have assumed the fields to be 
Gaussian in this calculation. See [13, HI] for discussions 
on the importance of non-Gaussian errors. 

For a CMB experiment, the fluctuations in the tem- 
perature field that is actually observed have contribu- 
tions from the primary CMB anisotropy, as well as from 
various secondary anisotropics and foregrounds such as 
the thermal SZ effect, the kinetic SZ effect, radio point 
sources, infrared point sources, diffuse Galactic fore- 
grounds etc, in addition to the instrumental noise [59j . 
Since the thermal SZ effect has a different frequency- 
dependence from the other components, by combining 
the data from multiple frequencies, the thermal SZ sig- 
nal can be extracted [see e.g. l60j. 

Here, for definiteness, we adopt the minimum- variance 
subtraction technique presented in (60j to extract the SZ 
signal, and, for simplicity, we include only the instrumen- 
tal noise. Therefore, the observed SZ power spectrum is 
given by 
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1^, 



(19) 



where we use the subscript "c" to label the quantities for 
different frequency channels; w — ((Tt6'fwhm/2cmb)~^, 
with ut the rms of the instrumental noise per pixel, and 
^'fwhm the full width, half-maximum (FWHM) of the 
beam, which is assumed to be Gaussian; is the Fourier 
transform of the beam profile, and Bj = exp[— £(£ -|- 
1)^fwhm/(^ In 2)]. Note we aim at extracting the SZ sig- 
nal in the RJ limit, so the observation at each frequency 
has been rescaled by the ratio of frequency dependence: 
s = —f{x)/2. Specifications of the CMB experiments 
considered in this paper are given in Table ID Since the 
instrumental noise is not correlated with the distribution 
of clusters, we have 



a 



yh 



C: 



■yh 



(20) 



Finally, we include a shot noise in the observed cluster 
power spectrum: 



^£,(a6,a'h') 



'~'i,(ab,a'b') 



1 



(21) 



Likewise, we can derive the covariance matrix for 
the stacked y-profile, (y) (^ab)i^i) ^ which is estimated by 
Fourier transforming the stacked SZ power spectrum and 
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TABLE I: CMB Experimental Specifications 



TABLE IL Optical Cluster Surveys 



Experiment ^^(deg^) 


u{GRz) 


^fwhm(') 




SPT[62l 2500 


95 


1.6 


26.3 




150 


1.1 


16.4 




220 


1.0 


85 


Planck[63l 41253 


30 


33 


5.5 




44 


24 


7.4 




70 


14 


12.8 




100 


10 


6.8 




143 


7.1 


6.0 




217 


5.0 


13.1 




353 


5.0 


40.1 




545 


5.0 


400.6 




857 


5.0 


18257.5 



Survey r2s(deg ) z-range 



SDSS 
DES 



7500 
5000 



0.1-0.3 
0.1-1.0 



LSST 


20000 


0.1-1.4 




III: Combined 


SZ and 


Optical 


Cluster 


Survey 


(deg 


f ■ 

) ■t'min 


^max 


SPT-DES 


2500 


7 ~ 


10000 


Planck-SDSS 


7500 


4 - 


- 3000 


Planck-LSST 


20000 


3 - 


- 3000 



then averaging over an annular area around 9i [see e.g. 
IgiI . for a similar derivation of the covariance for the cor- 
relation functions of the lensing field]. The covariance 
between (y)(a6)(^') ^^'^ ^v) (a'b')^^j) '^^ found to be 



tdl 



Mm) 



X Jo(^e,) (crC'^^aM'V) + C^U^fl'^')) ' (22) 



where Jo{i9i) is the average of Jo{(d) over the annulus 
around 9i. However, due to the exponential form of the 
noise term for C™, Eq. ([22]) leads to divergent results. 
To resolve this apparent problem, we instead calculate 
the covariance for the stacked beam-smoothed SZ pro- 
file, which is actually a direct observable. To calculate 
the beam-smoothed y-profile, we only need to multiply 
the integrand of Eq. ([T7)) by to account for the beam 
convolution, and to get the covariance, we only need to 
multiply that of Eq. ((22|) by _B|, which removes the ex- 
ponential divergence. Note that the covariance given by 
Eq. does not vanish when 9i ^ 6*^, i.e. the y-profiles 
at different angular bins are correlated with one another 
even when the field is Gaussian. 



the survey parameters. In this paper, we consider the 
case that an optical survey, overlapped with a CMB sur- 
vey, is used to construct a catalog of clusters. Such joint 
CMB and optical surveys that we are interested in in- 
clude, e.g., DES and SPT, Subaru HSC and ACT, and 
SDSS or LSST and Planck. 

As for the CMB surveys, specifications of each CMB 
experiment considered in this paper are given in Table ID 
Table |ll] gives the survey parameters for the optical sur- 
veys. Although an optical survey alone suffers from in- 
completeness of cluster find ing , which may lead to a com- 
plicated selection function 12 , |64 1 , we here assume that 
the optical survey can find all the clusters in the spec- 
ified ranges of mass and redshift. Obviously this is too 
optimistic, and one should keep in mind that the follow- 
ing results give the best available case for the stacked SZ 
measurements. For a more realistic study, the selection 
function for a given survey needs to be properly taken 
into account. Table Hill gives the parameters for the joint 
CMB and optical surveys, including the overlapping area, 
the minimum and maximum multipoles probed. 



B. Forecasts 



III. FORECASTS OF THE STACKED SZ 
MEASUREMENTS 

In this section, we estimate and discuss the detectabil- 
ity of the cross-correlation between the SZ map and the 
cluster distribution from ongoing and upcoming com- 
bined CMB and optical surveys. 

A. Survey Parameters 

To forecast the dctcctability of the stacked SZ signals 
from ongoing and upcoming surveys, we need to specify 



In Figures [T] and m we show the stacked SZ power spec- 
tra, scaled to the RJ limit, and the stacked i/-profiles for 
clusters in different redshift and mass bins. The dashed 
and dot-dashed curves in each panel are the one- and 
two-halo term contributions, respectively, while the solid 
curve represents the total. As can be seen, in Fourier 
space, the one-halo term dominates on small scales, while 
the two-halo term becomes more important on large 
scales, especially for the bins of clusters with lower masses 
and lower redshifts. Similar results hold in real space, ex- 
cept that the one-halo term is clearly truncated at about 
three times the virial radius of the most massive halo in 
each bin, as our gas profile is truncated at r = 3rvir (see 
§ lira)) . Note that we have converted the angular sepa- 
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FIG. 1: The cross power spectra of the SZ signal in the RJ band and the distribution of galaxy clusters in different redshift and 
mass bins. From left to right, the three columns are for clusters in the redshift bins of [0.2,0.3], [0.7,0.8] and [1.4,1.5], and from 
top to bottom, the rows are for clusters in the mass bins of [lO" *, 10" °], [lO"'^, lO^'* '*], [lO" '', 10" ®]/i"^Mq, respectively. 
The dashed and dot-dashed lines are the one- and two-halo term contributions, while the solid lines show the total. The boxes 
around each curve show ilcr statistical uncertainties in the power spectrum measurement expected for a 2500 deg^ SPT-DES 
like survey, calculated as we describe in Sec. IIIDI 



ration 9 to the projected distance separation rj^ by us- 
ing the (physical) angular diameter distance at the mean 
redshift of the bin. Thus, for a thin redshift slice, the 
stacking method probes the average projected pressure 
profile for the clusters stacked. In a statistically isotropic 
universe, the average 3D pressure profile is considered to 
be spherically symmetric, and therefore one can convert 
the projected pressure profile to obtain the 30 profile by 
using the Abel's integral [see e.g. Eq. (1) in[55|. 

The self-similar model predicts that the amplitude of 
the Compton y-parameter profile yo scales vifith the clus- 
ter's virial density pvir and virial mass Afvir as j/o 
Pvir-^vir (scc the Appendix). Since pvir increases with 
z, ?/o increases with z as well, in addition to its increas- 
ing with Mvir. These trends agree with what we see in 
Figures [T] and [2l the one-halo term for the cross cor- 
relation gets greater for cluster bins with higher z and 
larger M, though the gas model we use is not exactly 



self-similar (broken by the dependence of the concentra- 
tion parameter c on A/vir and z). As the clusters become 
more massive and their rcdshifts get smaller, the angles 
they extend become larger, hence in real space, the one- 
halo term terminates at a larger angular scale, while in 
Fourier space, it peaks at a smaller i. 

The error boxes around each curve in Figures [1] and 
[2] show the \a statistical uncertainties for measuring the 
stacked SZ signals at each multipole (radial) bin. These 
are expected for a 2500 deg^ SPT-DES like survey. Note 
the error boxes at different radial bins for the stacked y- 
profile are correlated. These figures show that the stack- 
ing method can lead to a significant detection of the SZ 
signal over a wide range of multipoles or radii. In par- 
ticular, it allows for exploration of the SZ signal at large 
radius, even around or beyond Tvir, as can be seen from 
Figure [21 We find the variance of Cf (Eq. [111) is mostly 
dominated by the term Cf^C'f'" . At £ ^ 100, it comes 
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0.1 1 10 0.1 1 10 0.1 1 10 

r,(Mpc/h) 



FIG. 2: Similar to Figure \T\ but here showing the stacked y-profiles for galaxy clusters in different redshift and mass bins, 
where the projected radius r± is converted from the angular separation 6 by using the angular diameter distance for the mean 
redshift of the cluster bin. The error boxes are for the beam-smoothed stacked profiles for the SPT-DES like survey, where we 
use the beam for the 150GHz SPT channel. Note that the errors at different radii are correlated. 



mainly from Cf^ /n^Jb ^ ^"^^ even smaller scales such as 
I > 10"^, instrumental noise dominates the error budget. 

To be more quantitative, we calculate the cumulative 
signal-to-noise square (S/N)'^ for measuring the stacked 
SZ signals for the clusters binned by their redshifts and 
masses, which in Fourier space is an integration over a 
given range of multipoles, defined as 




(23) 

The (S/N)'^ is equivalent to the information content 
available from a Fisher matrix analysis [65i l66j : given a 
template for the interested power spectrum, it gives the 
significance that the amplitude of the power spectrum is 
detected to be non-zero. We notice that the (S/N)'^ is 
proportional to the survey area ilg (assuming a simple 
survey geometry). 

In Figures [3][5l we plot the contours of the S/N for 
the SPT-DES, Planck-SDSS, and Planck-LSST like sur- 



veys, respectively. Again note that we have assumed 
all the clusters in a given mass and redshift bin can be 
found from the survey, i.e. 100% completeness and efh- 
ciency of cluster finding. As we noticed before, the am- 
plitude of the cross correlation increases when the clus- 
ters in the bin have larger masses or higher redshifts. 
At the same time, the number of clusters in the bin de- 
creases with their masses, and also decreases with their 
redshifts, though it increases with the redshift at first. 
While the S/N increases when the signal is higher, it 
decreases when the shot noise is larger. The result is 
that the S/N peaks for some cluster bin at interme- 
diate redshift and intermediate mass ~ the sweet spot, 
as can be seen from these three figures. From SPT to 
Planck, the resolution of the instruments gets worse, and 
the sweet spot shifts to lower redshift and higher mass, 
where the clusters extend larger angles such that their 
profiles can be better resolved. In particular, it would be 
interesting to notice that the arcminute-resolution CMB 
experiment such as SPT has the sweet spot at around 
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z ~ 0.6 and A/ ~ 2 x W^h-^MQ. From SDSS to LSST, 
the SjN increases due to the larger survey area. Note 
SjN oc ^ /sky , provided the probed i ranges are similar. 

Recently, the Planck team has shown their high signif- 
icance detection of the SZ signals for the SDSS maxBCG 
clusters [2^, down to the low mass systems. Their anal- 
ysis is different from what we are doing in this paper. 
Instead of stacking the SZ maps around the center of 
each cluster, they derived the integrated SZ signal Y500 
for each cluster by employing a multi-frequency matched 
filter, and then calculated the weighted mean of the sig- 
nal. The SjN of their detection ranges roughly 4^15 
for the cluster bins that are relatively narrower in mass 
and wider in redshift than ours. Despite the differences 
in analysis, their results are less optimal than ours. This 
may be because: the Planck results are from its first-year 
survey, so the instrumental noises are larger than their 
expected values, which are what we adopt in this paper; 
our results neglect the impurity and incompleteness of 
the maxBCG cluster catalog, hence tend to be optimal. 
Interestingly, we notice the SjN of the Planck's results 
also peaks at some intermediate mass ~ 2 x 10^'' ^."^Mq 
(corresponding to A^200 ^ 30), which is smaller than what 
we find by roughly a factor of 2. 




FIG. 4: Similar to the previous plot, but for a Planck-SDSS 
like survey. Note here the contours are displayed out to much 
higher redshift regions than those available from the SDSS, 
which is z:[0.1,0.3]. 
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FIG. 3: Contours of the total signal-to-noise ratio (S/N) ex- 
pected for measuring the cross-correlation between the SZ 
effect and the distribution of galaxy clusters assuming a SPT- 
DES hke survey (see Eq.[23]for the definition of the S/N). The 
clusters are binned in redshift with = 0.1 and in mass with 
Alog(Mobs) = 0.2, respectively. 



C. Residual Noise Contamination 

Our estimation for the noise in the extracted SZ sig- 
nal so far includes only the instrumental noise. This is 
optimistic, since there are various sources of contamina- 
tion to the SZ effect, e.g. the point sources (both the 
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FIG. 5: Similar to FigureEl but for a Planck-LSST like survey 

Poisson distributed and clustered components) |67| - 
[69j . and the primary CMB fiuctuations, which may leave 
some residual components even after the SZ signal has 
been maximally extracted. In this subsection, we con- 
sider two types of residual noise: those from the point 
sources and from the primary CMB, and see how the 
SjN of the stacked SZ signal is degraded. 

In the following, wc give a simple estimate on the im- 
pact of these two types of residual noise on the measure- 
ment of the stacked SZ signal. The primary CMB fluctu- 
ations do not correlate with the distribution of clusters. 
Therefore, we assume the residual CMB does not con- 
taminate the stacked SZ signal, but adds to its noise. We 
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also assume the residual CMB lias the same power spec- 
trum as the primary CMB, for simplicity. As for the con- 
tamination from point sources, we need to consider two 
populations: radio sources such as galaxies hosting active 
galactic nuclei (AGNs), and infrared sources such as the 
dusty star-forming galaxies. Although radio sources tend 
to reside in clusters, their contamination to the SZ signal 
is unlikely to be significant (10% at most for low-z clus- 
ters) [zO] ■ On the other hand, infrared sources dominate 
the point source contamination to the SZ power, with 
both a Poisson-distributed and a clustered component. 
However, it is likely that the majority of contamination 
comes from galaxies at high redshifts such as 2; ^ 2 [69j . 
which do not correlate with the massive halos at z ^ 1.5 
as are considered in this paper. Hence we assume that 
the residuals from point sources do not contaminate Cf ^, 
but add to C^^ . Specifically, we assume the residuals are 
Poisson distributed, and estimate their impact by simply 
doubling the instrumental noise. 

The results of our calculation are shown in Figure [51 
These are the S/N contours for a SPT-DES like survey. 
In the upper right panel, we include the residual noise 
from the point sources, in the bottom left, we include 
that from the primary CMB, while in the bottom right, 
we include both. To easily sec the effects of these resid- 
ual noises, we show in the upper left panel our results 
before, i.e., with instrumental noise only. As is expected, 
the S/N is reduced when either of the residual noises is 
included. According to our calculation for the residual 
power spectra, the primary CMB causes a larger reduc- 
tion in the S/N than the point sources, which perhaps 
indicates the necessity to remove the CMB. From Fig- 
ure [1] we can also see, on scales where the primary CMB 
dominates (£ <J 3000), there is non- negligible contribu- 
tion to the total S/N. Finally, we find, even in the most 
pessimistic case we study, the S/N for the detection of 
the stacked SZ signals for the less massive clusters fortu- 
nately remains substantial ( <; 20). 

IV. DISCUSSION 

In this paper, we have developed the halo model based 
method to compute the stacked SZ profile of the galaxy 
clusters, measured by cross-correlating the SZ map and 
the clusters' distribution, and the method to estimate 
the expected signal-to-noise {S/N) ratio for its detec- 
tion from a given survey. Such stacked SZ signals can 
be obtained, e.g., by combining the arcminute-resolution 
CMB experiments with optical surveys such as the com- 
bination of the DES and SPT, or the Subaru HSC and 
ACT surveys. 

Although the measured CMB signal is by nature two- 
dimensional, cross-correlating the CMB map with clus- 
ters of known redshifts allows for a tomographic recon- 
struction of the SZ signals as a function of redshift. In 
addition, simultaneous reconstruction of the SZ signals as 
a function of cluster mass is also allowed, if the clusters' 
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FIG. 6: Similar to Figure [S] but using different models for 
the noise in the SZ field. The upper-left panel is the same 
as Figure O shown here for comparison, which includes only 
the instrumental noise, the upper-right has the residual noise 
from point sources added in, the lower-left has that from the 
primary CMB added in, while the lower right has both (see 
text for the details). 



masses are also known. This would be very useful for 
carrying out the precision cosmology with cluster based 
experiments as well as studying the cluster astrophysics, 
in particular, its time evolution [see, e.g. .[2^. 

The stacked SZ signals will be very complementary 
to the stacked lensing signals [l^, [s^ , if both methods 
study the same catalog of clusters. One of the crucial 
sources of systematic in cluster cosmology is from the 
mass-observable relation. To attain the precision of clus- 
ter cosmology, we need a well-calibrated mass-observable 
relation. Combining the stacked SZ and lensing signals 
for the same catalog of clusters will provide a power- 
ful means of studying the mean relation between cluster 
mass and SZ signal as a function of redshift and mass 
(more precisely, the optical richness) from the observa- 
tion. As we have shown, the stacking method can sig- 
nificantly reduce the statistical errors, and allows one to 
probe the SZ or lensing signals for the low mass clus- 
ters or to large angular scales [also see[l3l- However, the 
scatters around the mean mass-observable relation still 
need to be understood, and detailed studies of the target 
cluster sample combining various observables will also be 
needed [e.g., seelTJ, for such a multi- wavelength study]. 

Combining the stacked SZ and lensing signals will be 
powerful for studies of cluster physics as a function of 
cluster mass and redshift. Due to the statistical sym- 
metry of the universe, the stacked pressure and mass 
profiles can be considered to be spherically symmetric 
in an average sense. Hence it would be of great interest 
to construct the cluster's baryon fraction profile obser- 
vationally, and then address at which radii the cosmic 
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mean baryon fraction is reached. Since the intra-cluster 
gas may be blown off to outer radii by various feedback ef- 
fects from e.g. AGN or supernova, which involve compli- 
cated astrophysical processes, such a model-independent 
method of studying the baryon fraction would be very 
useful [see [7l|, [t^, for the related discussion] . 

The SZ signals from galaxy clusters would be smaller 
if there is a significant fraction of non-thermal pressure 
support from random gas motions, whose existence could 
help to explain the apparent conflict between the ampli- 
tude of the observed SZ power spectrum with theoretical 
prediction (20| . Such non-thermal pressure support is ex- 
pected to be more significant at the outskirts for given 
mass-scale halos [2l| . Again the stacking method is very 
powerful for studying the SZ signals at such outer radii, 
which otherwise would be very difficult due to the faint 
signals. 

Throughout this paper, we have been using the rela- 
tively simple analytical model from for the cluster 
pressure profile, given as a function of halo mass and 
rcdshift. It is known this model overestimates the mea- 
sured SZ power spectrum amplitude (by about a factor 
of 2), which however would not alter our conclusion of a 
significant S/N for the stacked SZ measurement that is 



evident from Figure 3-5. However, a more careful study 
will be definitely needed, e.g ., by using more realistic 
pressure profiles [see, e.g.,|73|, or hydrodynamic simula- 
tions, in order to realize the genuine power of the stacking 
method. 
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APPENDIX 



In the self-similar model, the gas profiles of all the clusters are considered identical with appropriate choice of 
normalization factors that are the powers of the global quantities e.g. mass Mvir, radius Tvir- For example, the 
temperature profile can be expressed as T(r)/Tvir oc /T(?'/''vir), with Tvir c>c Mvir /''vir (from the virial theorem) and 
/^(r/rvir) a universal function. Similarly, the pressure profile can be expressed as P(r)/Pvir oc fp{r /r^ir), with 
Pvir oc PvirTvir and fp{r/r^i^) a universal function. Hence, the y-profile is given by 

(X PvirMvir/y(6'/6'vir), (24) 

where 0vir ~ 'Vir/c^A, with cIa the angular diameter distance, and fy is a universal function. As can be seen, the 
amplitude of the self-similar y-profile is proportional to pvir-^^vir- 



